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Abstract 

Let L be a Schrodinger operator (^V — A{x)) 2 + V{x) with periodic 

CO ■ magnetic and electric potentials A, V, a Maxwell operator V x -4-rVx 

>, £\ x ) 

in a periodic medium, or an arbitrary self-adjoint elliptic linear partial 

,_i- ' differential operator in M n with coefficients periodic with respect to 

C^) ■ a lattice T. Let also S be a finite part of its spectrum separated by 

J , gaps from the rest of the spectrum. We consider the old question of 

existence of a finite set of exponentially decaying Wannier functions 

qq ■ Wj(x) such that their T-shifts Wj n {x) = Wj(x — 7) for 7 € T span 

(^ . the whole spectral subspace corresponding to S in some "nice" man- 

£> I ner. It is known that a topological obstruction might exist to finding 

exponentially decaying Wj^ that form an orthonormal basis of the 
spectral subspace. This obstruction has the form of non-triviality of 
certain finite dimensional (with the dimension of the fiber equal to the 
number m of spectral bands in S) analytic vector bundle A5 over the 
n-dimensional torus. It was shown by G. Nenciu in 1983 that in the 
presence of time reversal symmetry (which implies absence of mag- 
netic fields), and if S is a single band, the bundle is trivial and thus 
the desired Wannier functions do exist. In 2007, G. Panati proved that 
in dimensions n < 3, even if S consists of several spectral bands, the 
time reversal symmetry removes the obstruction as well, if one uses 
the so called composite Wannier functions. It has not been known 
what could be achieved when the bundle is non-trivial (which can 
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occur for instance in presence of magnetic fields or for Chern insula- 
tors). Let r be the type of the bundle Ag, i.e. the number of open 
sub-domains over which it is trivial (for the trivial bundle r = 1, and 
r never exceeds 2 n , where n is the dimension of the coordinate space) . 
We show that it is always possible to find a finite number / < rm 
(and thus m < I < 2 n m) of exponentially decaying composite Wan- 
nier functions Wj such that their T-shifts form a 1-tight frame in the 
spectral subspace. Here 1-tight frame is a redundant analog of an 
orthogonal basis, which appears in many applications, e.g. to signal 
processing and communication theory. This appears to be the best 
one can do when the topological obstruction is present. The number 
/ is the smallest dimension of a trivial bundle containing an equiva- 
lent copy of As- In particular, / = m if and only if A$ is trivial, in 
which case an orthonormal basis of exponentially decaying composite 
Wannier functions is known to exist. 



1 Introduction 

Wannier functions, along with Bloch waves, play an important role in solid 
state physics [1], photonic crystal theory (see, e.g., [21,26] for general discus- 
sion of photonic crystals and [7,34,36] for Wannier function applications), 
and in fact in any study that involves periodic linear elliptic equations. They 
are heavily involved in description of electronic properties of solids, theory 
of polarization, photonic crystals, numerical analysis using tight-binding ap- 
proximation, etc. (see, e.g. [1,5,7,23,28-32,34,36,38,39] and references 
therein). As it is said in [29], strongly localized Wannier functions "are the 
solid-state equivalent of "localized molecular orbitals"..., and thus provide 
an insightful picture of the nature of chemical bonding." It is in particular 
important to have the Wannier functions decaying as fast as possible. That 
is why the question of possibility of choosing a finite number of families of 
exponentially decaying Wannier functions that form an orthonormal basis in 
the spectral subspace corresponding to an isolated part S of the spectrum has 
been intensively considered in physics literature, at least since the paper by 
W. Kohn [23], who showed that this is possible in ID. However, the problem 
becomes non-trivial in 2D and higher dimensions. And indeed a topological 
obstruction arises there which is not present in ID. Namely, existence of 
such a family is known [35] to be equivalent to triviality of certain analytic 
vector fiber bundle that we will call A^ (see (19)). It was shown by G. Nenciu 



[31] in 1983 (see also [18]) that in the presence of time reversal symmetry, 
and if S is a single spectral band, the bundle is trivial and thus the desired 
family of Wannier functions does exist. In 2007, G. Panati [32] proved that 
in dimensions n < 3, even if S consists of several spectral bands, the time 
reversal symmetry removes the obstruction as well. In this case, though, one 
should naturally resort to the so called composite, or generalized Wannier 
functions that correspond to a finite family of bands, rather than to a single 
band. The activity in this direction seems to be high and still increasing 
[29,39]. 

It has not been known what could be achieved when the bundle is non- 
trivial (which can happen in presence of magnetic fields or for Chern in- 
sulators). In this case, a family of Wannier functions with the described 
properties cannot exist. One can avoid the topological obstruction relaxing 
the exponential decay to the square integrability of the Wannier functions, 
but this amounts to a decay that is considered to be way too slow. There 
is not much one can improve on this decay without striking the topological 
obstruction (e.g., a slow decay such as summability of L 2 norms over the 
shifted copies of the Wigner-Seitz cell is already impossible, see Theorem 
5.4). 

In this text, we show existence of what appears to be the best exponen- 
tially decaying Wannier family one can expect in the case of a non-trivial 
bundle Ag. In order to formulate the main result, we need to introduce some 
notions. 

Definition 1.1 

• A frame in the Hilbert space H is a set of vectors ipj such that for any 
vector f G H the following inequality holds: 



^ii/n 2 <Ei(^)i 2 <5 



The frame is said to be A-tight ; if B = A in the inequality above. In 
particular, for a 1-tight frame one has 



£l(/^. 



X 



The reader should notice that there is no condition of linear independence 
on the vectors forming a frame, so frames are usually redundant (and thus 
are not bases). 



It is known (see, for instance [17]) that 1-tight frames are exactly orthog- 
onal projections of orthonormal bases from a larger Hilbert space and play a 
substitute of such bases for redundant systems of vectors. If a 1-tight frame 
is linearly independent, it in fact forms an orthogonal basis. 

Frames in general, and tight frames in particular, have become recently 
a standard tool in signal processing and communication theory, in particular 
due to their numerical stability and robustness with respect to noise (e.g., 
[3,4,8,9,11,12,22] and references therein). 

We also need one more notion: 

Definition 1.2 (e.g., [20]) Let A be a (locally trivial) vector bundle over 
a connected base T . 

• The bundle is said to be of a finite type, if there is a finite open 
covering {£^}i<j<t of the base T such that the bundle A is trivial over 
each of Uj 

• The type r of the bundle is the minimal number of elements of such a 
covering. 

• We will call the dimension (also called the rank) of the vector bundle 
A the dimension of its fibers. 

Remark 1.3 

• Vector bundles over a compact base (e.g., over a torus) always have a 
finite type. 

• The type is equal to 1 if and only if the bundle is trivial. 

• The type of any vector bundle over an n-dimensional torus does not 
exceed 2 n . 

We can now formulate the following main result of the paper: 
Theorem 5.7 Let L be a self-adjoint elliptic V -periodic operator in IR ra , n > 
1 and S C K be the union of m spectral bands of L. We assume that S is 
separated from the rest of the spectrum. Let also t be the type of the obstacle 
bundle As over the n-dimensional torus. Then there exists a finite number 
I < rm (and thus m < I < 2 n m) of exponentially decaying composite Wan- 
nier functions Wj(x) such that their shifts Wj n := Wj(x — 7), 7 G T form a 



1-tight frame in the spectral subspace H$ of the operator L. This means that 
for any f(x) G H$, the equality holds 



\f(x)\ dx = 2_^ | / f(x)w j:J (x)dx\ . (1) 

3,1 R n 

The number I G [m, 2 n m] is the smallest dimension of a trivial bundle con- 
taining an equivalent copy of As- In particular, I = m if and only if As is 
trivial, in which case an orthonormal basis of exponentially decaying compos- 
ite Wannier functions exists. 

The number / of the Wannier functions Wj has to exceed the number of 
bands, unless there is no topological obstacle in the form of non-triviality of 
As- This creates a redundant system of functions, and thus the orthonormal 
property is not achievable. The 1-tight frame property is the best analog of 
orthonormality one can get in this case. For instance, it allows the control 
of the L 2 norms in terms of the projections onto the Wannier system. 

The structure of the paper is as follows: Section 2 surveys the main 
notions of the Floquet-Bloch theory that will be needed for the rest of the 
text (they can be found in most solid state texts or in [24-26, 33]); in Section 3 
we introduce the Wannier functions and discuss their properties of interest; 
Section 4 contains an overview of the known results concerning the single 
band case; Section 5 moves to the composite bands and composite Wannier 
functions. It also contains the proof of the main result of the paper. The first 
half of the paper provides a self-contained survey of known results, which is 
done for the benefit of the reader. This part also sets up all preliminary 
results and notions that are needed for the proof of the main result of the 
text. 



2 Main notions and preliminary results 

Let L(x, D) be a bounded from below 1 self-adjoint elliptic operator in M n . 
The specific nature of the operator will be irrelevant (e.g., matrix operators, 
such as Dirac or Maxwell can be allowed). One can thus consider, without loss 
of generality, one's favorite periodic operator, e.g. the Schrodinger operator 
(iV — A(x)) 2 + V(x) with real periodic magnetic and electric potentials 



This restriction is not essential and can be removed. 



A, V. Some conditions need to be imposed on the potentials to define a self- 
adjoint operator L in L 2 (R n ) (e.g., [10,33]). However, for what follows these 
details are insignificant, and thus one can safely assume "sufficiently nice" 
coefficients of L. 

Let T be a lattice in M. n , i.e. the set of integer linear combinations of 
vectors of a basis a\,...,a n in M. n . (Due to the general form of L that we 
allow, no generality will be lost if the reader assumes that T is just the integer 
lattice I/ 1 .) The coefficients of the operator L are assumed to be periodic 
with respect to the shifts by vectors 7 of T. We will be using a fundamental 
domain W (e.g., the Wigner-Seitz cell in physics literature) of the lattice I, 
i.e. such a domain that its T-shifts cover the whole space with only boundary 
overlap. 

We denote by T* the reciprocal (or dual) lattice to T. It lives in the 
dual space to M ra , but if an inner product (•, •) in R n is fixed, then T* can 
be realized in the same space M n as consisting of all vectors k, such that 
(k,7) G 2ttZ for all 7 G T (e.g., if T = Z n , then T* = 2vrZ n ). A fundamental 
domain B for T will be fixed (e.g., the first Brillouin zone in physics). The 
quotient M n /r* forms a torus that we will denote T*. It can also be realized 
as B with the boundary points identified according to the action of T*, and 
thus folded into a torus. In particular, T*-periodic functions on W" 1 are more 
naturally identified with functions on T*, rather than with ones on B. 

The periodicity of the spectral problem Lu = Xu with respect to the lat- 
tice group T suggests to use the group Fourier transform. This immediately 
leads to the well known under various names both in physics and mathe- 
matics [1,24-26,33] transform, which we will call here the Bloch-Floquet 
transform: 

/(*) -> /(*, x) := J2 /(* + 7)e-* 7 - (2) 

Here k is a real (we will also need complex values of k) n-dimensional vector, 
which we will call quasi-momentum (its other common names are crystal 
momentum and Bloch momentum). Assuming that / decays sufficiently fast 
(e.g., if / has bounded support, or is in L 2 (R n )), there is no convergence 
problem 2 . It is straightforward to check that for any (even complex) quasi- 
momentum k, the function f(k, x) is T*-periodic with respect to k and of the 



2 In fact, like the Fourier transform, the Bloch-Fourier transform can be extended in a 
distributional sense to a much wider class of functions, e.g. [25]. 



Bloch (also called Floquet) form with respect to x, which means that 

f(k,x) = e ik -*v k (x), (3) 

where Vk(x) is T-periodic. Thus, the values x G W and k G B are sufficient 
for determining the whole function f(k,x). In particular, one can consider 
f(k,x) as a function f(k,-) on B with values in a space of functions on 
W. Due to r*-periodicity, it is more natural (and important for us for what 
follows) to consider f(k, •) as a function on the torus T* rather than on 
B. Formally, this means the change of variables from the quasi- momenta 
k — (ki, . . . , k n ) to Floquet multipliers z = e lk := (e lk ' ai , . . . , e lk ' an ). Then 
the torus T*, which is the image of the space of all real quasi- momenta, 
gets imbedded into the complex space C n as the unit torus {z||zj| — l,j — 
1, . . . , n}. The space of all complex quasi-momenta becomes the set (C\0) n C 
C n that contains all complex vectors with non-zero components. We will also 
need some complex neighborhoods of the space of real quasi-momenta and 
of the torus T*, defined for a given a > 0: 

V a = {k G C n | |Im k-dj\<a,j = l,...,n}, (4) 

and its image under the transform k i— > z 

Q a = {z = (z 1 ,...,z n ) G Q\e~ a < \zj\ < e a , j = l,...,n}. (5) 

Here {ctj} is the basis of T mentioned before. 

If / G L 2 (R n ), the series (2) converges in the space L 2 (B,L 2 (W)) of 
square integrable functions on B with values in L 2 (W). 

Another simple remark which we will need is the following: T-shifts of / 
result in multiplication by exponents in the Bloch-Floquet images. Namely, 
if oj G T and f u (x) = f(x — a;), then 

Uk,x) = e- ik -"f{k,x). (6) 

In the case of constant coefficient linear differential operators, plane waves 
are generalized eigenf unctions, and thus the Fourier transform diagonalizes 
the operators. In the periodic case, Bloch functions (3) and Bloch-Floquet 
transform (2) play the roles of plane waves and Fourier transform correspond- 
ingly. 

In order to formulate simple but crucial for us analogs for Bloch-Floquet 
transform of the standard Plancherel, inversion, and Paley- Wiener theorems 
for the Fourier transform, we need to introduce some notions first: 



Definition 2.1 The space L 2 (W l ) consists of all functions f G L 2 oc (M. n ) such 
that for any < b < a the following expression is finite: 

Mf) '■= SU P ll/l|L2(w+ 7 )e 6171 < oo. (7) 

7er 

This space is equipped with the natural topology defined by the semi-norms 

So, this is the space of functions which are locally squarely integrable, and 
whose local L 2 norm decays exponentially when the domain is shifted to 
infinity. When we will be discussing the exponential decay of a Wannier 
function w, this will always mean exponential decay of ||w||£2( W+7 ) with 
respect to 7. 3 

We will also denote by A(Q,H), where Q is an n-dimensional complex 
domain and H is a Hilbert space, the space of all if -valued analytic functions 
on Q, equipped with the topology of uniform convergence on compacta. 

We are now able to formulate the promised basic properties of the Bloch- 
Floquet transform (e.g., [25]): 

Theorem 2.2 

1. If f G L 2 (M. n ) and K C M. n is a compact, then the series (2) con- 
verges in the space L 2 (T* , L 2 (K)) . Moreover, the following equality 
(Plancherel theorem) holds: 

* m ;ll/(Mll!«, w « = /ll/MllL (w ,^. <«) 

B T* 

where d'k is the normalized to total measure 1 Lebesgue measure on B ; 
and d ' z is the normalized Haar measure on T* . 

2. For any a G (0, 00], Bloch-Floquet transform 
is a topological isomorphism of the space L 2 ,(W a ) onto A(Q a ,L 2 (W)). 



3 If the coefficients of the operator L are good enough, one might improve the norm 
in which the decay is observed. For instance, if L maps functions that locally belong to 
the Sobolev space iJ™/ 2 + e (with some e > 0) into functions that are locally in U n / 2 - r + e j 
where r is the order of the operator (usually r = 2 in all physics applications), one can 
get pointwise decay. 



3. For any f G L 2 (R") the following inversion formula holds: 

f(x)= I f(k,x)d'k, xeW 1 . (9) 



The first statement claims that the Bloch-Floquet transform is an isome- 
try between the natural Hilbert spaces, the second shows that exponential 
decay transforms into analyticity in a neighborhood of the torus T* (a Paley- 
Wiener type theorem), and the third one provides inversion formulas for the 
transform. 

Due to periodicity of the eigenfunction equation Lu = Aw, the Floquet 
transform block diagonalizes it. This leads to the well known (see [1,24,25, 
33] and references there) description of the band-gap spectral structure of 
the operator L. Namely, let L(k) be the operator L acting on the Bloch 
functions (3) with a fixed quasi-momentum k. If we denote (for real k) the 
eigenvalues of L in nondecreasing order as 

Ai(ife) < \ 2 (k) < . . . i-> oo, (10) 

then 

1. band functions Xj(k) are continuous, r*-periodic, and piece-wise an- 
alytic in k; 

2. if Ij is the finite segment that is the range of Xj(k) (the jth band), 
then the spectrum of L is 



a(L) = \jL 



j- 
:i 



Analytic properties of the band functions and the corresponding Bloch 
eigenfunctions <f)j(k, x) (i.e., L<pj = \j<pj) are considered in detail in [25, 41]. 
In particular, one can always choose <pj(k,-) as a piece-wise analytic L 2 - 
function on B with values in L 2 (W), whose norm in L 2 (W) is almost ev- 
erywhere constant and can be chosen equal to 1. We will also assume below 
that <f)(k, •) is r*-periodic with respect to k. 

For constant coefficient operators, the basis of delta functions is dual 
under Fourier transform to the basis of plane waves and is also very useful. An 
attempt of inventing in this vein an appropriate for periodic problems analog 
of delta functions immediately leads to the so called Wannier functions. 

9 



3 Wannier functions 

Let <f>j(k,x) be a Bloch eigenfunction (not necessarily normalized) corre- 
sponding to the band function Xj(k) (the choice of such an eigenfunction 
is non-unique, even when the eigenvalue is simple). We will assume that 
<f>j G L 2 (T*, L 2 (W)), which is known to be always possible ([25,41] and ref- 
erences there). 

Definition 3.1 The Wannier function Wj(x) corresponding to the Bloch eigen- 
f unction (f)j(k,x) is 

Wj(x)= </>j(k,x)d% xeW 1 . (11) 



Comparing this definition with (9), one concludes that the Wannier func- 
tion Wj is just the inverse Bloch-Floquet transform of <f)j(k,x), and corre- 
spondingly <f)j(k,x) is the Bloch-Floquet transform of wf 

(f>i(k, x) = w](k, x) = J2 w j( x + 7)e~^ 7 - (12) 

The analog of this property for the usual Fourier analysis is that plane waves 
can be obtained as Fourier transforms of delta functions. 

Equality (12) enables one to rephrase the definition of a Wannier function 
as follows: 

w(x) is a Wannier function for a periodic operator L, if its Bloch-Floquet 
transform w(k, x) for any k E M" is an eigenfunction of the corresponding 
Floquet operator L(k). 

It is often useful to have a Wannier function Wj(x) with mutually orthog- 
onal lattice shifts Wj(x — 7). In fact, these shifts are often considered as 
different Wannier functions. It is immediate to check that 

Wj (x + -f) = f e ik ^(j) j (k,x)d'k for 7 e T. (13) 

B 

One also wants to have normalized Wannier functions. In these directions, 
simple answers are provided by the following corollary of (6) and Theorem 
2.2: 

Corollary 3.2 

10 



1. The Wannier function Wj(x) belongs to L 2 (M. n ) and 



Wj (x)\ 2 dx= / U 3 {k,-)\\ 2 LHw) d'k. (14) 



2. Functions Wj^(x) := Wj(x — , ~f) are mutually orthogonal for 7 G T if and 
only if the eigenfunction <pj(k,x) in (17) has a k-independent norm in 
L 2 (W). 

The first claim follows from the first statement of Theorem 2.2. The second 
one follows from (13) and the fact that the only periodic functions on B that 
are orthogonal to all exponents e* fc ' 7 with 7 e r, 7 7^ 0, are constants. 

If the family <f)j(k, x) is not normalized, then the shifts of the correspond- 
ing Wannier function are not mutually orthogonal. However, if the norm 
||(/>j(/c, -)||l2( W ) never vanishes, one can normalize it to 



to get a new Wannier function 



</>j(k 



3Vi 



Uj(k, -)IU 2 (w) 



w'Ax) = I tpj(k,x)d'k. 



that already has orthogonal shifts. 

The most important property of Wannier functions in many problems of 
physics is their decay. The already established square integrability is far from 
being sufficient, and exponential decay is often desired. 

Considering the Bloch eigenfunction <f>(k, x) as a vector- valued function 
of the quasi- momentum k, it is clear (as a direct consequence of Theorem 
2.2) that smoothness of <pj(k,x) with respect to k translates into decay of 
Wj(x). (Since periodicity of <f>j(k) is assumed, smoothness also means match- 
ing the values and all derivatives across the boundaries of Brillouin zones.) 
In particular, the following simple statement holds (see [25, Section 2.2] and 
Theorem 2.2 above): 

Lemma 3.3 

1. If Y2 \\ w j\\l 2 cw+7) < °°; then <f)j(k,-) is a continuous L 2 (W)-valued 
function on T*. 

11 



2. Infinite differentiability of<ftj(k,-) as a function on T* with values in 
L 2 (W) is equivalent to the decay of ||wj||x,2( W+7 ) faster than any power 

o/ItI f° r l e r - 

3. Analyticity of <pj(k,-) as a function on T* with values in L 2 (W) is 
equivalent to the exponential decay of \\wj\\L^(w+'y)- 

Theorem 2.2 provides more details about the relation between the decay and 
analyticity domain in the last statement of the Lemma. 

Thus, it is important to address the issue of choosing <f)j(k, •) as smooth 
with respect to k as possible. We will look at this problem first for a single 
spectral band. 

4 A single band case 

Dealing with a single band Xj(k) and the eigenfunction branch <f)j(k, •), one 
encounters (at least in the mult i- dimensional case, unlike ID [23]) a problem 
whenever band functions cross, since then smooth continuation of the eigen- 
function and the eigenvalue through the crossing might not be (and usually 
is not) possible [42]. Thus, one should assume a band that does not intersect 
other bands. In this case, although the band function A.,- (A;) itself is periodic 
and analytic with respect to k in a neighborhood of the real space, choosing 
an analytic branch of the eigenfunction that never vanishes might still be 
impossible. In fact, even continuity of this function might not be achiev- 
able. Indeed, finding an analytic and periodic <pj(k,-) such that the norm 
\\(f)j(k, -)IU 2 (w) never vanishes, means finding a non-vanishing periodic ana- 
lytic section of the one-dimensional periodic analytic complex vector bundle 
over M. n : 

A r .= \J ker(L(k) - Xj(k)I) (15) 

Here ker A means the kernel of the operator A, i.e. the space of solutions 
Au = 0. 

Taking into account r*-periodicity and considering Floquet multipliers 
z = e tk instead of quasimomenta, one can reinterpret this bundle as one over 
the torus T*. The latter view is the one we will adopt here, while keeping 
the same notation for the bundle. 

Due to one-dimensionality of the bundle, existence of a continuous non- 
vanishing section over T* of Aj is equivalent to its triviality. The bundle, 

12 



however, might be non-trivial, which will be a topological obstruction to 
existence of not only analytic, but even continuous non-vanishing section. 

One might wonder whether existence of an analytic section (and thus 
of an exponentially decaying Wannier function with orthogonal shifts) faces 
more obstructions than just the existence of a continuous one (which would 
lead to a much slower decaying Wannier function with orthogonal shifts). As 
the following proposition states, this happens to be not the case 4 . Before 
formulating it, let us notice first that if a > is small enough, the band 
function Xj(k) extends analytically into the neighborhood T> a of the real 
space M. n C C n (or to its periodized version Q a , as we will implicitly assume) 
without colliding with other band functions. Then the bundle Aj naturally 
extends analytically to Q a : 

K ha := |J ker(L(A;) - \,(k)I), z = e lk . (16) 

z£il a 

Proposition 4.1 The bundle Aj over T* is topologically trivial if and only 
if for a small a > the bundle Aj yQ over Q a is analytically trivial. 

Proof. We notice first of all that Q a is a product of one-dimensional complex 
domains, and thus is a Stein manifold (e.g., [16]). Then a Grauert theorem 
(an incarnation of "Oka's principle") applies that says that topological and 
analytic triviality of A_j j0 are equivalent [13-15]. Since T* is a deformation 
retract of Q a , the triviality of Aj >a over Q a is equivalent to triviality of Aj 
over T*. □ 

Thus, one needs to worry about topological obstructions only. There is 
one important case when these obstructions do not materialize. This is when 
the operator L has real coefficients. Then the complex conjugate u to any 
solution u of the equation Lu = Xu with real A is also a solution. In particu- 
lar, if u(x) = e lk ' x v(x) is a Floquet solution with the quasimomentum k, then 
u(x) = e~ lk ' x v(x) is a Floquet solution with the quasimomentum —k. Thus, 
in the Floquet structure of L the symmetry k h- > —k is present. In physics 
this corresponds to the time reversal symmetry and holds for Schrodinger 
operators with real electric potentials, but breaks down in the presence of a 
magnetic field (since it contributes imaginary first order terms to the opera- 
tor). 

The following result is proven in [31]: 



4 This statement has been well known and is implicitly present in [31,35]. 
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Theorem 4.2 [31] If the coefficients of the self-adjoint elliptic V -periodic 
operator L are real and Xj{k) is an analytic band function without crossings 
(for real k) with other band functions, then the bundle Aj j0 (defined for small 
a > 0) is analytically trivial. Hence, there exists a corresponding exponen- 
tially decaying Wannier function Wj(x) whose Y -shifts Wj n {x) = Wj(x — 7) 
are mutually orthogonal and normalized in L 2 (M. n ). 

If the band function Xj does not have crossings with other band func- 
tions, in dimensions n > 1 the corresponding band Ij of the spectrum still 
can overlap with other bands. Assume now that this does not happen and 
the band Ij = [dj,bj] is surrounded by gaps separating it from the rest of 
the spectrum. Then there is a decomposition of the space L 2 (R") into the 
orthogonal sum of L-invariant subspaces 

L 2 {R n ) = H i ®HJ-, 

such that the spectrum of L in Hj coincides with the band Ij, while the 
spectrum of L H ± is the rest of the spectrum of L. In other words, Hj is the 
spectral subspace of L corresponding to Ij. Then Theorem 4.2 leads to 

Corollary 4.3 [31] In the presence of time reversal symmetry, if the spectral 
band Ij does not intersect other bands of the spectrum of L, then there ex- 
ists an exponentially decaying Wannier function Wj(x) such that its T -shifts 
Wj n (x) = Wj(x — 7) form an orthonormal basis in Hj. 

Proof. According to Theorem 4.2, there exists an analytic, r*-periodic in k, 
and normalized function (f>j(k, •) G Ker(L(k) — Xj(k)I). According to Corol- 
lary 3.2, its inverse Bloch-Floquet transform gives a Wannier function Wj with 
normalized and orthogonal T-shifts Wj r/ . It only remains to show that these 
functions are complete in the spectral subspace Hj to which they all belong. 
Indeed, after Bloch-Floquet transform, the subspace Hj becomes the space 
of all vector functions of the form f(k)<f>(k, •), where f(k) is any r*-periodic 
scalar function such that J \\f(k)\\ 2 d'k < 00. On the other hand, finite lin- 

B 

ear combinations of Wannier functions J^ a^Wj n transform into functions 

r 

Y^ « 7 e* fc ' 7 I (f)(k,-). Since, according to the standard Fourier series theory 
r / 

[37], trigonometric polynomials I Yl o 7 e lfe ' 7 J are r*-periodic and dense in 
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L 2 (B), this proves completeness of the set of Wannier functions Wj n in Hj. 

D 



5 Multiple bands and composite Wannier func- 
tions 

Suppose now that a set S consists of m > 1 bands of the spectrum (with 
overlaps and crossings allowed) and is separated from the rest of the spectrum 
by spectral gaps. An example would be the part of the spectrum from its 
bottom till the first gap. Let us denote the union of these bands by S (a 
composite band). A natural question to ask now is the following: are there 
analogs of Theorem 4.2 and Corollary 4.3 for this case? Again, there exists an 
orthogonal decomposition into L-invariant subspaces H$ ® Hg, where H $ is 
the spectral subspace of L corresponding to the union of bands S. After the 
Floquet transform, functions / G H$ will correspond to functions f(k, •) that 
for each fc£B belong to the spectral subspace Hs,k £ L 2 (W) of the Floquet 
operator L(k) that corresponds to the set S. 

Since the band functions that correspond to the bands in S might cross, 
the attempts to find analytic non-degenerate families of eigenf unctions 4>(k, x) 
are expected to be futile in general. However, one can relax the requirement 
that <f)(k,x) is an eigenf unction of L(k) and request only that <p(k, •) G H$ t k 
for each real k. This leads to the following notion: 

Definition 5.1 A composite (or generalized) Wannier function w(x) 
corresponding to an isolated from the rest of the spectrum union S of bands 
is the one that can be represented as 

w{x) = <p{k,x)d'k, (17) 

T* 

where <p(k, •) G H$ t k for all real k. 

Let us assume that S consists of m bands. One can ask the question: do 
there exist m exponentially decaying generalized Wannier functions Wj such 
that their shifts Wj n (x) = Wj(x — 7), 7 G T form an orthonormal basis in the 
spectral subspace Hs of the operator L ? 

In order to attempt an answer, let us consider a simple contour C in the 
complex plane that surrounds the set S and separates it from the rest of the 
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spectrum of L. It does not intersect the spectrum of L(k) for any real k, and 
thus, by a perturbation argument, does not intersect it for all quasimomenta 
k with sufficiently small imaginary parts. Hence, if one considers the Tri- 
dimensional spectral projector Ps{k) for L(k) in L 2 (W) that corresponds to 
this contour 

Ps(k) = ^-<f((I-L(k))- 1 d(, (18) 

2m J 

c 

the projector will be analytic with respect to A; in a (r*-periodic) neighbor- 
hood of the space of real quasimomenta k. Thus, its range forms an analytic 
r*-periodic vector bundle 

As = [jPs(k)(L 2 (W)) (19) 

over a neighborhood T> a of the real space. As before, taking into account 
periodicity and considering Floquet multipliers z = e tk instead of quasimo- 
menta, one can reinterpret this bundle as one over a neighborhood Q a of 
the unit torus T*. The latter view is the one we will adopt here, without 
changing notations for the bundle. 

Taking into account our previous discussion and the requirement that the 
number of the "mother" Wannier functions Wj is assumed to be equal to the 
number m of bands in S, one obtains the following 

Lemma 5.2 [32] Let S be the union of m bands that is isolated from the 
rest of the spectrum. A family of m exponentially decaying generalized Wan- 
nier functions Wj such that their shifts Wj yl (x) = Wj(x — 7), 7 G T form an 
orthonormal basis in the spectral subspace Hg of the operator L exists if and 
only if the bundle Ag over Q a is analytically trivial. Equivalently (due to 
Oka 's principle), this happens if and only if the restriction of A5 to T* is 
topologically trivial. 

The following complement to Theorem 4.2 requires non-trivial topological 
arguments and was proven in [32]: 

Theorem 5.3 [32] Let the coefficients of the self-adjoint elliptic V -periodic 
operator L in M. n , n < 3 be real (and thus the time reversal symmetry holds) 
and S C M. be the union of m spectral bands of L. We assume that S is 
separated by gaps from the rest of the spectrum. Then the bundle A5 is 
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analytically trivial. Hence, there exist m exponentially decaying generalized 
Wannier functions Wj(x) whose T -shifts Wj n (x) = Wj(x — 7) form an or- 
thonormal basis in H$. 

The nice Theorems 4.2 and 5.3 give the desired positive answer in physical 
dimensions n < 3 and any m and for any dimension n when m = 1 , whenever 
the time reversal symmetry is present. However, without this symmetry (or 
in dimensions 4 and higher, when m > 1) one cannot guarantee triviality 
of the bundle and thus existence of the required generalized Wannier basis 
in the spectral subspace S. One can wonder what can be achieved without 
triviality of the bundle. According to Lemma 5.2, if the bundle A5 is non- 
trivial, one should sacrifice some of the properties guaranteed in Theorem 
5.3: exponential decay of generalized Wannier functions, their number (i.e., 
one might need more Wannier functions than the number m of the bands), 
orthogonal basis property, or several of these. One might try to weaken the 
exponential decay condition to avoid redundancy in the system of Wannier 
functions and thus to keep the basis property. The next simple statement 
shows that in the presence of the topological obstacle, one cannot go far in 
this direction, since even the weak decay condition (20) below is already as 
strong as the exponential one. 

Theorem 5.4 Let L be a self-adjoint elliptic V -periodic operator in M. n ,n > 
1, and S C H. be the union of m spectral bands of L. We assume that S is 
separated from the rest of the spectrum. Suppose that there exists a family of 
m composite Wannier functions Wj(x) whose T-shifts form an orthonormal 
basis in Hs and such that for each j = 1, . . . , m the following sum is finite 

^ ||wj||l2(w+ 7 ) < 00. (20) 

7er 

Then the bundle A5 is trivial and thus there exists a system of m exponen- 
tially decaying Wannier functions Vj that have the same basis property. In 
particular, power decay as 

iKiiL 2 (w +7 )<c(i + i7ir n - e 

for some e > (or a stronger requirement that \w(x)\ < C(l + \x\)~ n ~ € ) 
already implies absence of the topological obstacle. 
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Proof. Indeed, according to Lemma 3.3, condition (20) implies continuity 
with respect to k of the corresponding Bloch functions <pj(k, •). Together with 
orthonormality condition, this gives a basis of continuous sections, and thus 
triviality of Ag. This, in turn, implies, as we have already seen, existence of 
the corresponding family of exponentially decaying Wannier functions. □ 

So, we would like to preserve the exponential decay. The next statement 
shows that one can do so at the cost of increasing the number of Wannier 
functions and thus sacrificing their orthogonality, while still having a com- 
plete system of functions. This is a preliminary simple result, since Theorem 
5.7 proved later improves significantly on the completeness statement and 
provides an estimate of redundancy of the system of Wannier functions. 

Theorem 5.5 Let L be a self-adjoint elliptic V -periodic operator in M. n ,n > 
1, and S C K. be the union of m spectral bands of L. We assume that S is 
separated from the rest of the spectrum. Then there exists a finite number 
I > m of exponentially decaying composite Wannier functions Wj(x) such 
that their shifts Wj n , 7 G T form a complete system in the spectral subspace 
Hs (i.e., the set of all their finite linear combinations is dense in the space). 

Proof. We will prove this statement in an equivalent formulation: there 
exists a finite number of analytic sections <pj(k, •), j = 1, . . . , / of the bundle 
As over Q a for some a > such that the sections <pj,k{k, •) := e lk ' 1 (f)j(k, •) for 
j — 1, ...,/, 7 G T form a complete set of functions in the space of L 2 -sections 
of the bundle Ag. It is a simple exercise to show that it is sufficient for that 
to have a finite family of analytic sections <pj such that they span the whole 
fiber of the bundle at each point z G T*. Recall now that the bundle As 
is trivial locally, i.e. in a neighborhood of each point. Let us fix two finite 
trivializing covers Uj and Vj <<= Uj , j = 1 , . . . , r of T* by complex domains 
such that their union is a Stein domain (e.g., the union can be made equal 
to fi a for some a > 0). Then on each Uj, by assumption, there is an analytic 
basis if)j tt (k),t = 1, . . . , m in the restriction of the bundle As onto Uj. Due to 
the Stein property, these holomorphic sections can be approximated with any 
precision uniformly on Vj by global (i.e., defined on the whole VL a ) analytic 
sections <f>j,t(k) (see, e.g. [19, Ch. VII] or [25, Theorem 1.5.9]). Doing this 
for each j, with a sufficiently close approximation, one achieves the desired 
property for the functions (f)jj- d 

In this result we seemingly have lost too much: we have an excessive 
(in comparison with the number of spectral bands involved) and uncontrol- 
lable number of composite Wannier functions, and besides the orthonormal 



property is lost, being replaced by a much weaker (and not very useful) com- 
pleteness. Well, one cannot avoid having more than m generalized Wannier 
functions. Indeed, the following statement holds: 

Proposition 5.6 If in the previous theorem either the number of Wannier 
functions is exactly equal to m, or the shifts of the Wannier functions Wj 
form an orthonormal basis, then the bundle As is trivial. 

Proof. Indeed, assume that the number of the functions is exactly m. 

It is not hard to notice that in this case the Floquet transforms <f)j of the 
Wannier functions must be analytic and at each point must span the whole 
fiber of the bundle. Since the fibers have dimension m, we conclude that <pj 
must form an analytic basis of the bundle, and thus the bundle is trivial. 

Suppose now that the functions Wj n are orthonormal. Normalization 
property, as we already know, means that the functions <pj(k, •) have a con- 
stant norm in L 2 (W). Checking orthogonality of Wj lt o to Wj 2a for ji ^ j 2 
and arbitrary 7 G T, one obtains that 

3 ik -' y ( ( f) jl (k, •), <f> j2 (k, -))v(yr)dk = 

for all 7 G T. This means that the whole Fourier series with respect to k of the 
function (4>j 1 (k, •), (f>j 2 (k, -))l 2 (w) onT* vanishes. Thus, (<pj 1 (k, •), (f>j 2 (k, -))l 2 (w) 
is equal to zero almost everywhere, and hence the functions <pj 1 , <pj 2 are or- 
thogonal for each k. Since the dimension of the fiber is m, we conclude that 
the number of the Wannier functions does not exceed m, and thus, according 
to the previous case, the bundle is trivial. □ 

So, if one cannot help with excessive number of Wannier functions, can 
one improve on the completeness statement? It seems that the answer must 
be negative, since an overdetermined system of vectors cannot be orthogonal. 
However, we have described in Definition 1.1 the notion of the so called tight 
frame of vectors (e.g., [17]) that replaces orthonormality in the redundant 
case. This is what allows us now to improve on the statement of Theorem 
5.5 and prove our main result. 

Theorem 5.7 Let L be a self-adjoint elliptic V -periodic operator in M n ,n > 
1 and S C K. be the union of m spectral bands of L. We assume that S 
is separated from the rest of the spectrum. Let also t be the type of the 
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obstacle bundle As over the n- dimensional torus T* . Then there exists a finite 
number I < rm (and thus m < I < 2 n m) of exponentially decaying composite 
Wannier functions Wj(x) such that their shifts Wj n := Wj(x — 7), 7 G T form 
a 1-tight frame in the spectral subspace Hs of the operator L. This means 
that for any f(x) G Hs, the equality holds 

\f(x)\ 2 dx = J2\ [ f(x)w~(x)dx\ 2 . (21) 

3,1 R n 

The number I G [m, 2 n m] is the smallest dimension of a trivial bundle con- 
taining an equivalent copy of As- In particular, I = m if and only if As is 
trivial, in which case an orthonormal basis of exponentially decaying compos- 
ite Wannier functions exists. 

Proof. Let us establish first of all the following auxiliary statement: 

Lemma 5.8 There exists a finite dimensional analytic sub-bundle A' of the 
trivial bundle Q a x L 2 (W) ; such that 

1. The fibers at the same point k of As and A' intersect only at zero. 

2. The analytic sub-bundle $ := A5 A' in Q a x L 2 (W) is trivial and 
has dimension I, with m < I < 2 n m as described in Theorem 5.7. 

There are several ways to show that this lemma holds. For instance, one 
way is to choose an abstract finite-dimensional analytic bundle A' such that 
added to A5 restricted to some Qf, (where Q is the closure of Q) it leads to 
a trivial bundle of dimension /. This is possible in the topological category 
[2,20] and thus, according to Grauert's theorem [15] and Stein property of 
Q a , also in analytic category. Using Gaussian maps, one can make sure (e.g., 
proposition 3.5.8 in [20]) that the dimension I of the direct sum As A' does 
not exceed rm, where m is the dimension of the fibers of the bundle A^ and 
r is its type (see Definition 1.2). In the case of an n-dimensional torus as the 
base, as it is mentioned in Remark 1.3, r does not exceed 2 n . Thus, one can 
always achieve / < 2 n m. 

Now let us consider a complementary to A5 infinite dimensional sub- 
bundle JF in Q a x L 2 (W), which always exists [42, Theorem 3.11] (and in our 
case can be constructed using the spectral projectors, as in (18)). The claim 
is now that there is an analytic sub-bundle in T that is isomorphic to A'. 
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Indeed, taking the direct sum of A' and of the trivial bundle with an infinite 
dimensional separable Hilbert space G as a fiber, one gets a Hilbert analytic 
fiber-bundle over a Stein manifold. According to Kuiper's theorem [2,27], 
this bundle is topologically trivial. Bungart's theorem [6] (an infinite dimen- 
sional version of Grauert's theorem [15] on fiber bundles on Stein spaces) 
implies that this sum is also analytically trivial. Since, again by the Kuiper 
and Bungart theorems, all analytic infinite dimensional Hilbert bundles with 
separable fibers over Stein bases are equivalent, we conclude that the direct 
sum is equivalent to J 7 , which embeds a copy of A' into T . This proves the 
Lemma. 

Let us now prove the statement of the theorem, using the sub-bundle 
A' as a tool. Consider the trivial /-dimensional bundle $ := As@A' C 
Q a x L 2 (W). Due to its triviality, one can find a basis of its analytic sections 
ipj(k),j = 1, . . . , /. What we need to do is to make this system orthonormal 
for real values of k (or equivalently for z G T*) without losing the analyticity. 
This is possible in our case. Indeed, one would like to use the standard 
Gram-Schmidt orthogonalization method, but in a version that preserves 
analyticity with respect to the quasi-momentum k (equivalently, with respect 
to the multiplier z). Consider the domain T> a , which is a neighborhood 
in C n of the space of real quasimomenta M n . It is clearly symmetric with 
respect to the complex conjugation k = (hi, . . . , k n ) i— > k — (hi, . . . , k n ). 
Let <f>(k) be an analytic function in T> a with values in a complex Hilbert 
space H (H = L 2 (W) will be our case) with the Hermitian scalar product 
(•, •) (which is thus anti-linear with respect to the second factor). Then the 
function (<f>(k) , <f>(k)) is analytic in V a and coincides with ||0(/c)|| 2 for real 
k. If now ||0(/c)|| does not vanish for real k, then one can normalize 4>(k) as 
follows: 

^(fc) ;= ffl (22) 



Then ip(k) is analytic in a neighborhood of the real space and has unit norm 
for real values of k. Notice that the expression in (22) can lose its analyticity 
at some point, due to the denominator, and thus the neighborhood T> a might 
need to shrink somewhat. An analogous analytic procedure can be applied 
to inner products by computing (cf>(k),ip(k)). Now this allows one to apply 
the Gram-Schmidt orthonormalization for real k in a manner that preserves 
analyticity in a complex neighborhood of W 1 . The same procedure applies in 
a neighborhood Q a of the torus T*. We can thus produce another analytic 
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basis family in $, which is ortho normal for real values of k. We can assume 
now that ipj{k) is already like that. Applying the spectral projector Ps{k) 
to this family of sections, we obtain an analytic family of sections (j)j(k) of 
the bundle As, which in each fiber form a 1-tight frame (as the orthogonal 
projection of an orthonormal basis in a larger space). Consider the two 
families of generalized Wannier functions: Vj that correspond to ipj and Wj 
that correspond to <f)j. The functions Vj n do not necessarily belong to the 
spectral subspace H$, but according to Theorem 2.2 and Corollary 3.2, they 
form an orthonormal basis in a Hilbert subspace H' that is larger than Hs- 
(This is the subspace that after the Bloch-Floquet transform produces the 
space of all L 2 -sections of the bundle $ over the torus.) Thus, the family 
Wj t y, as the orthogonal projection of Vj n from H' onto Hg, form a tight frame 
in H s . □ 

6 Conclusion 

1. We have shown that for any isolated finite part S of the spectrum 
of an elliptic self-adjoint operator periodic with respect to a lattice T, 
there exists a family of I exponentially decaying composite (generalized) 
Wannier functions such that their T-shifts form a tight frame in the 
spectral subspace H s that corresponds to S. The number / is equal to 
the smallest possible dimension of the fiber of a trivial vector bundle 
containing the obstacle bundle A$. In particular, it does not exceed 
Tin, where r is the type of the bundle As and m is the number of 
spectral bands constituting S. Thus, m < I < 2 n m. The numbers I 
and m are equal, and thus there is no redundancy in the constructed 
set of Wannier functions, if and only if As is trivial, and then the shifts 
of these Wannier functions form an orthonormal basis in Hs- 

2. It was also shown that an attempt to relax the decay condition in 
presence of the topological obstruction does not work well, since even 
a very slow decay condition (20) hits the same topological obstruction. 

3. The results apply to a general class of elliptic self-adjoint periodic oper- 
ators, including systems such as Dirac operators and Maxwell operators 
in periodic media. Magnetic translations can also be treated the same 
way. 
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4. The results also apply in an abstract situation of a periodic elliptic 
self-adjoint operator on an abelian covering of a compact manifold, 
graph, or quantum graph. Neither formulations, nor proofs require any 
modifications in this case. 
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